A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11(11.09. 1;12. 05. 1;11.01.2;11.06.1;11.09.4;02.12.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 



o 
o 
o 

(N 

X> 
<D 
IX, 



> 
o 

(N 
(N 
O 

o 
o 



Oh' 
I 1 

O 1 



X 



A vestige low metallicity gas shell surrounding the radio 
galaxy 0943-242 at z = 2.92 

L. Binette 1 , J. D. Kurk 2 , M. Villar-Martm 3 , H. J. A. Rottgering 2 

1 Instituto de Astronomi'a, UNAM, Ap. 70-264, 04510 Mexico, DF, Mexico (e-mail: binette@astroscu.unam.mx) 

2 Leiden Observatory, P. O. Box 9513, 2300 RA, Leiden, The Netherlands 

3 Department of Physical Sicences, University of Hertfordshire, College Lane, Hatfield Herts, ALIO 9AB, England 
Received / Accepted 



Abstract. Observations are presented showing the dou- 
blet CivAA1548, 1551 absorption lines superimposed on 
the C IV emission in the radio galaxy 0943-242. Within 
the errors, the redshift of the absorption system that has 
a column density of Nciv = 10 14 5±01 cm -2 coincides 
with that of the deep Lya absorption trough observed by 
Rottgering et al. (1995). The gas seen in absorption has 
a resolved spatial extent of at least 13kpc (the size of the 
extended emission line region). We first model the absorp- 
tion and emission gas as co-spatial components with the 
same metallicity and degree of excitation. Using the infor- 
mation provided by the emission and absorption line ratios 
of C iv and Lya, we find that the observed quantities are 
incompatible with photoionization or collisional ionization 
of cloudlets with uniform properties. We therefore reject 
the possibility that the absorption and emission phases are 
co-spatial and favour the explanation that the absorption 
gas has low metallicity and is located further away from 
the host galaxy (than the emission line gas). The larger 
size considered for the outer halo makes plausible the pro- 
posed metallicity drop relative to the inner emission gas. 
In absence of confining pressure comparable to that of the 
emission gas, the outer halo of 0943-242 is considered to 
have a very low density allowing the metagalactic ionizing 
radiation to keep it higly ionized. In other radio galaxies 
where the jet has pressurized the outer halo, the same gas 
would be seen in emission (since the emissivity scales as 
n H ) and not in absorption as a result of the lower filling 
factor of the denser condensations. This would explain the 
anticorrelation found by Ojik et al. (1997) between Lya 
emission sizes (or radio jet sizes) and the observation (or 
not) of H I in absorption. The estimated low metallicity for 
the absorption gas in 0943-242 (Z ~ O.OIZq) and its pro- 
posed location -outer halo outside the radio cocoon- sug- 
gest that its existence preceeds the observed AGN phase 
and is a vestige of the initial starburst at the onset of 
formation of the parent galaxy. 
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1. Introduction 

Very high redshift (z > 2) radio galaxies (hereafter 
HZRG) show emission lines of varying degree of ex- 
citation. In virtually all objects, the Lya line is the 
strongest and is usually accompanied by high excitation 
lines of C iv AA1549, C in] A1909, He n A1640 and, at times, 
NvA1240 (Rottgering et al. 1997 and references therein). 
An important characteristic of the emision gas is its spa- 
tial scale. The sizes of the Lya emission region range from 
15 to 120 kpc (van Ojik et al. 1997). 

Most ground work on HZRG is performed at rather 
low resolution (~ 20A) to maximize the probability of 
line detection and the S/N. However a very potent discov- 
ery was made by van Ojik et al. (1997, hereafter v097) at 
much higher resolution, that of extended Hi absorption 
gas. In effect, out of 18 HZRG spectra taken at the un- 
usually high resolution of ~ 1.5-3A, v097 found -in 60% 
of the objects- deep absorption troughs superimposed on 
the Lya emission profiles. Furthermore, out of the 10 ra- 
dio galaxies smaller than 50 kpc, strong Hi absorption is 
found in 9 of them. The absorption gas appears to have a 
covering factor near unity over very large scales, namely 
as large as the underlying emission gas. 

The current paper addresses the problem of the ion- 
ization state of both the absorption and the emission gas 
as well as the interconnection between the two. The main 
justifications behind this work are the following: HZRG 
are probably the progenitors of the massive central clus- 
ter galaxies (Pentericci et al. 1999) and as such are an 
important means by which we can study large ellipticals 
and their environment at such high redshift, a time not 
so long after, or even during their formation. Further- 
more, the extended gas as detected in Civ (see below) 
is chemically enriched and therefore represents debris of 
past intense stellar formation periods and is interesting 
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to study in their own right. What is the fate of such gas? 
How quickly has the enrichment of this large scale gas pro- 
ceeded? Will this gas be heated up into a hot wind and 
enrich the intergalactic X-ray gas in cluster of galaxies? 
Will it on the contrary condense into sheets or conden- 
sations? A better understanding of the various gaseous 
phases which co-exist in high redshift objects would help 
anwering these questions. 

To determine the physical conditions of the absorption 
gas, new observations were carried out at the wavelength 
of Civ and Hen in 0943-242, the first radio galaxy re- 
ported to show large scale absorption troughs (Rottgering 
et al. 1995, hereafter R095). The new spectrum shows 
the Civ absorption doublet at the same redshiftQ Zd as 
the Lya absorption trough (R095). Clearly and surpris- 
ingly the gas in absorption is highly ionized and probably 
of comparable excitation to the gas seen in emission. 

This paper is structured as follows. We first present 
observations which show Civ in absorption in 0943-242 
(Sect. ||). In Sect. || we derive a ratio (r) relating the 
observed emission and absorption quantities which de- 
pends somewhat on the ionization fraction of H but not 
explicitely on the C/H metallicity ratio. At first, we pos- 
tulate that the emission and absorption gas components 
are co-spatial and share the same excitation mechanism 
and physical conditions and proceed to model T with a 
one-zone equilibrium photoionization model. We improve 
on the model using a stratified photoionized slab. As the 
observed ratio cannot be reproduced even in the case of 
collisional ionization, we discuss in Sect. 4.1 two alter- 
native interpretations of this significant discrepancy. We 
demonstrate the many advantages of the winning scenario 
in which the absorption gas is further out and of much 
lower density, pressure and metallicity than the emission 
gas. 



2. Observations of C iv (and Lya) in absorption 
in 0943 242 

2.1. Earlier observations of 0943-242 at z e = 2.92 

The low resolution spectrum of 0943-242 shown in R095 
and discussed also in van Ojik et al. (1996) displays the 
characteristic emission lines of a distant radio galaxy: 
strong Lya, weaker C IV, He II and possibly C in] . This ob- 
ject was also observed at intermediate resolution (1.5 A) by 
R095 in the region of Lya with the slit positioned along 
the radio axis. The initial discovery of extended absorp- 
tion troughs was based on this latter spectrum which we 
reproduce in Fig. [l]. 
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Fig. 1. An expanded plot of the Lya spectral re, 
tained by R095. The Hn emission gas redshift 
2.9233 ± 0.0003 and the main absorber of column 



$ion ob- 
is z e = 
Nhi = 
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19.0±0.2 



cm- 2 lies at z a = 2.9200 ± 0.0002. 
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We will distinguish between absorption and emission red- 
shifts using subscripts, as in z a and z e , respectively. 



Fig. 2. The full-resolution AAT spectrum showing the 
CivAA1548, 1551 and HeiiA1640 lines. 



2.2. New observations of C IV and Hell at 
intermediate resolution 

With the objective of providing constraints on the abun- 
dances and kinematics of the gas in 0943-242, sensitive 
high-resolution spectroscopic observations centered at the 
C iv and He n lines were performed at the Anglo Aus- 
tralian Telescope (AAT) on 1995 March 31 and April 1 
under photometric conditions and with a seeing which var- 
ied from 1" to 2". The RGO spectrograph was used with a 
1200 grooves mm -1 grating and a Tektronix 1024 2 thinned 
CCD, yielding projected pixel sizes of 0.79" x 0.6A. The 
projected slit width was 1.3", resulting in a resolution as 
measured from the copper-argon calibration spectrum of 
1.5A FWHM; the slit was oriented at a position angle of 
74°, i.e. along the radio axis (as in R095). 
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The total integration time of 25000s was split into 
2 x 2000s and 7 x 3000s exposures in order to facilitate re- 
moval of cosmic rays. Exposure times were chosen to en- 
sure that the background was dominated by shot noise 
from the sky rather than CCD readout noise. Between 
observations the telescope was moved, shifting the object 
slit by about 3 spatial pixels, so that for each exposure 
the spectrum was recorded on a different region of the 
detector. The individual spectra were flat-fielded and sky- 
subtracted in a standard way using the long-slit package 
in the NO AO reduction system IRAF. The precise offsets 
along the slit were determined using the position of the 
peak of the spatial profile of the C iv and He II lines. Using 
these offsets, the images were registered using linear inter- 
polation and summed to obtain the two-dimensional spec- 
trum. The resultant seeing in the final two-dimensional 
spectrum, measured from two stars on the slit, was 1.5" 
FWHM. The corresponding FWHM of C iv emission along 
the slit was 2.2", giving a deconvolved (Gaussian) width 
of 1.6" or 12kpc. Within the errors, this is the same as 
that found for Lya emission by R095. 

The two-dimensional spectrum was weighted summed 
over a 7 pixel (5") aperture to obtain a one-dimensional 
spectrum. In Fig. |2| we show the AAT data in the form of 
a full-resolution spectrum. 

2.3. Profile fitting of the emission and absorption Lya 
and C IV lines 

One deep trough is observed in the Lya emission line 
(Fig. 0) which was interpreted as a large scale Hi ab- 
sorber by R095. In addition there are a number of weaker 
troughs, presumably due to weak H I absorption. Fit- 
ting the emission line by a Gaussian and the H I ab- 
sorption by Voigt profiles, R095 infer a column den- 
sity Nhi of io 19 0±0 - 2 cm -2 for the deep trough, a red- 
shift z a — 2.9200 ± 0.0002 and a Doppler parameter b 
of 55 ± 5 kms" 1 . For the three shallow troughs, they 
find Nhi ranging from 10 13 8 to 10 141 cm -2 and b rang- 
ing from 7 to 100 kms" 1 . The redshift difference of the 
absorbers relative to systemic velocity when converted 
into inflow/outflow velocities indicate values not exceeding 
800 kms . Because at the bottom of the main trough no 
emission is observed, the covering factor of the absorbing 
gas must be equal or larger than unity over the complete 
area subtended by the Lya emission, indicating that the 
spatial scale of the absorber exceeds 13kpc. This work will 
concern only the deep absorption trough. 

To parameterize the C iv profile we have assumed that 
the underlying emission line is Gaussian, with Voigt pro- 
files due to the C IV doublet absorption superimposed. We 
used an iterative scheme that minimizes the sum of the 
squares of the difference between the model and the ob- 
served spectrum, thereby solving for the parameters of the 
model (e.g. Webb 1987, v097). Initial values were assumed 
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Fig. 3. An expanded plot of the full-resolution spectrum 
with the fit superimposed (solid line). 

Table 1. Parameters for the Gaussian and Voigt profile 
fits 



Emission C iv He II 

Offset (10" 11 erg cm" 2 s" 1 ) 0.29 ± 0.01 0.32 ± 0.07 

Peak (10" 17 erg cm" 2 s" 1 ) 1.90 ± 0.1 1.75 ± 0.2 

Position of peak (A) 6078.2 ± 0.5 6434.5 ± 0.5 

z e 2.9247 ± 0.0003 2.925 ± 0.001 

FWHM (kms" 1 ) 1430 ± 50 1025 ± 45 

Absorption C iv 

~ a 2.9202 ± .0002 

b (kms" 1 ) 45 ± 15 

Nciv (cm " 2 ) 10 i4.5±o.i 

Position of 1 st trough (A) 6068.2 

Position of 2 nd trough (A) 6078.3 



for the shape of the Gaussian profile and the redshift of 
the absorber. 

In Fig.|| we show a portion of the spectrum with the 
model fits superimposed. The Gaussian fitted to the C iv 
emission line peaks at z e = 2.9247 ± 0.0003 and has a 
FWHM of 29 ± 2A. We have corrected all wavelengths to 
the vacuum heliocentric system (~+1.13A) before com- 
puting the redshifts. The two troughs in this figure corre- 
spond to the CivAA1548, 1551 doublet produced by the 
same absorption system. Therefore, within the fitting pro- 
cedure, the wavelength separation and the ratio of the two 
profiles' depths are fixed by atomic physics while the two 
values for b are set to be equal. The fit gives for the lo- 
cation of the bottoms of the two troughs A = 6068.2 and 
6078.3A resulting in a redshift of 2.9202 ± 0.0002. Within 
the errors this redshift is equivalent to that of the main 
Hi absorber and in the subsequent analysis we will as- 
sume that the Lya and C iv absorption gas belongs to the 
same absorber. We derive a Doppler parameter b for the 
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doublet of 45 ± 15kms 1 and a column density Nciv of 
-^q14.5±o.i cm -2 as sumrnar i ze d in Table |l|. 

As expected, He n appears only in emission without 
any absorption since it is not a resonance line. Parame- 
ters for the He n emission profile were obtained by fitting 
a Gaussian using the same iterative scheme (see Fig. f 
in Rottgering & Miley 1997). The peak is positioned at 
z e = 2.925 ± 0.001 and has a FWHM of 22 ± 2k. The 
fitted parameters of the emission and absorption profiles 
are presented in Table [l]. We recall that the FWHM of the 
Lya emission profile is 1575 ± 75 kms -1 (v097), signifi- 
cantly larger than that of He n (see Table [l]) . Inspection 
of the various profiles in Fig. |l| and Fig. || (or Fig. ||) sug- 
gests the presence of an excess flux on the blue wings of 
all the emission profiles. Combining information from all 
the emission lines, our best estimate of the emission gas 
redshift is z e = 2.924 ± 0.002. 

2.4- Velocity shear and subcomponents 

To investigate whether there is any velocity shear in the 
C iv emission profile we fitted spatial Gaussian profiles to 
the emission line as a function of wavelength. In Fig. |^ we 
show the wavelength maxima of these spatial profiles and 
a line fitted through these points. The spatial profile of the 
C iv emission spectrum is displaced by 0.2", corresponding 
to a displacement of 1.5 kpc, over a wavelength range of 
50A. R095 measured a comparable shift for Lya of 1.8 
kpc| although it appears that the latter displacement is 
due to a far more pronounced and abrupt difference in 
locations of the Lya peak on both sides of the absorption 
trough. As Fig.|J shows, the peaks of C IV emission form a 
wavy line. We believe the velocity shear in the Civ profile 
to be less significant than the shear in the Lya profile. 
We cannot rule out that the small velocity shear might 
be masking a possible break up of the absorption regions 
into a few saturated absorption components of smaller b. 

A concern about the determination of Nciv is the pos- 
sibility that that there exist subcomponents in the absorp- 
tion systems that have high column densities but low b 
values and are, therefore, not acounted for whenever indi- 
vidual velocity subcomponents are not resolved. Although 
we cannot strictly exclude this possibility, we adopt the 
stand of Jenkins (1986) and Steidel (1990a) who, using 
extensive absorption line studies, argue that this is un- 
likely to be the case, at least for C iv, and that a single- 
component curve-of-growth analysis can be used to infer 
total columns although the inferred effective b value has 
no physival meaning in terms of temperature. It is inter- 
esting to note that the physical conditions inferred from 
the C iv fit are fully consistent with the observed ratio of 
the doublet (since both troughs are equally well fitted). 

2 This new value of 0.33 ± 0.06 pixels x 0.74 arcsec/pixel = 
0.2442 arcsec x 7.36 kpc/arcsec = 1.80 ± 0.33 kpc is to be pre- 
ferred to that quoted by R095 of 2.5 kpc. 
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Fig. 4. The relative location of the peak of the C IV emis- 
sion at constant wavelength as a function of wavelength. 
The line is a weighted fit to these peaks. The zero offset 
is arbitrary. 

If the underlying continuum was flat, the Nciv column 
and the b value we infer would imply a theoretical ratio of 
equivalent widths of Wo(1548)/W (1551) = 1.4, which is 
where the curve of growth just begins to leave the linear 
part (Steidel 1990a). Clearly the Njji column might be 
susceptible to a larger error since Lya is saturated. With 
these caveats in mind, we will assume in the following 
analysis that the adopted columns do not lie far off from 
reality. 

3. A simple model for the ionized gas in emission 
and absorption 

Our initial hypothesis is that the absorption gas is a sub- 
component of the emission gas, sharing the same excita- 
tion mechanism and metallicity. We discuss the physical 
conditions of such gas and proceed to calculate an observ- 
able quantity, T, against which to compare the information 
provided by the Lya and Civ lines in 0943-242. 

3.1. Relation between the ionized absorption and 
emission components 

The C iv and Lya lines are both resonant lines and there- 
fore prone to be seen in absorption against a strong un- 
derlying source. This property has consequences for the 
emission gas as well. In effect, for a geometry consisting 
of many condensations for which the cumulative covering 
factor approaches unity, the resonant line photons must 
scatter many times in between the condensations before 
they can escape. In this case, the emerging flux of any 
resonant line from a non uniform distribution of gas will 
not in general be an isotropic quantity but will depend on 
geometrical factors and on the relative orientation of the 
observer, a point which we now develop further. 
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We propose that some kind of asymmetry within the 
emission gas distribution can explain how a fraction of the 
ionized gas can be seen in absorption against other nearby 
components in emission. Let us suppose that the emission 
region is composed of low filling factor ionized gas conden- 
sations which are denser (therefore brighter) towards the 
nuclear ionizing source. In this picture, the Lya or C IV 
photons are generated within and escape from such con- 
densations, after which they start scattering on the surface 
of neighboring condensations until final escape from the 
galaxy (we assume that the cumulative covering factor is 
unity). Let us now suppose an asymmetry^ in the global 
distribution of the outer condensations respective to the 
plane of the sky. In this case, the total number of scatter- 
ings on neighboring condensations before final escape will 
differ depending on the perspective of the absorber. Since 
for an observer situated on the side with an excess of con- 
densations many of the resonant photons would have been 
'reflected' away, we expect that the reduced flux would ap- 
pear as an absorption line at the same velocity as that of 
the condensations responsible for reflecting away the res- 
onant photons. The outer condensations (responsible for 
the absorption) must necessarily be of lower density in or- 
der to be of negligible emissivity respective to the inner 
(denser and therefore brighter) emission gas, otherwise the 
outer gas would out-shine in emission! 

We should point out that for a density of the absorp- 
tion gas as high as 100 cm -3 as argued for in v097, such 
a gas cannot be photoionized by the metagalactic back- 
ground radiation which would be much too feeble to pro- 
duce C iv. The ionization to such a degree of the absorp- 
tion gas is in itself puzzling. We adopt as working hypoth- 
esis that it is -similarly to the emission gas- photoionized 
by the AGN or by the hard radiation from photoionizing 
shocks. 

Finally, the fact that both the absorption and emission 
gas contain a significant amount of C +3 argues in favor of 
a common geometry and excitation mechanism for the gas, 
the underlying hypothesis behind the calculations devel- 
oped below. 



3.2. The observable quantity T 

The quantities determined from observation of 0943-242 
are the following: the emission line ratio measured by 
Rottgering et al. (1997) is = 0.194. We adopt the 

value of 0.17 following estimation of the missing flux due to 
the absorption troughs. As for the absorption gas, the Hi 
and Civ column densities are 10 19 cm" 2 and 10 14 5 cm~ 2 , 
respectively, as discussed in Sect. ||. These four quantities 
carry information on the three ionization species H , H + 



3 The asymmetry would take place either in space or in ve- 
locity domain or both. 



and C +3 . We define the ratio T as the following product 
of the emission and absorption ratios: 
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where Nhi / Nciv is the ratio of the measured absorption 
columns. If, as postulated above, the gas responsible for 
absorption is simply a subset of the line emitting gas, the 
ratio r does not explicitly depend on the abundance of 
carbon as shown below. 

3.3. The simplest case of an homogeneous one-zone slab 

To compute T, in a first stage let us consider an homo- 
geneous slab of thickness L of uniform gas density, tem- 
perature and ionization state to represent both the gas 
in emission and in absorption. Ignoring any peculiar scat- 
tering effects, the emission line ratio j° IV is given by the 
ratio of the local emissivities jciv I jLya since the slab is 



homogeneous. For the emissivity of the C IV line, we have 



Anjciv — 8.63 10 6 hvc iv n e nciv 
x-^exp(-to OT /ff)/\f 



(2) 



(Ostcrbrock 1989) where T is the temperature, flciv the 
collision strength of the combined doublet, u>i the statisti- 
cal weight of the ground state and hvciv the mean energy 
of the C IV excited level. For the Lya emissivity, we have 

^j Lya = hv Lya n e n HII aJJ (T) 



(3) 



where cS^J is the effective recombination coefficient rate 
to level 2p of H (Osterbrock 1989). By putting the tem- 
perature dependence and all the atomic constants in the 
function f(T), the emission line ratio becomes: 



Iciv _ Z^ nl nHi]civ 



LLya 



nnVnii 



f(T) 



(4) 



where n# is the total hydrogen density, Z^T 11 the carbon 
abundance relative to H of the emission gas, rjciv the frac- 
tion of triply ionized C and ynii the ionization fraction 
of H. 

The ratio of column densities Nhi I Nciv can be writ- 
ten as: 



Nhi 



n H x H i 



Nciv Z c bs n H n C iv 



(5) 



where xhi is the neutral fraction of H inside our homo- 
geneous slab and Z c bs the carbon abundance of the ab- 
sorption gas. As we are testing the case which equates the 
absorption gas with the emission gas, then Z c bs — Z^V 11 . 
We denote as T the product of the two calculated ratios: 
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We note that T is not directly dependent on either the 
abundance of C or on its ionization state. It is, however, 
dependent on the temperature and on the ionization state 
of H through the ratio^ fj^fj- To compute this ratio, it is 
necessary to postulate an excitation mechanism. For this 
purpose, we have used the code MAPPINGS ic (Binette, 
Dopita & Tuohy 1985; Ferruit et al 1997) to compute la- 
under the assumption of either collisional ionization or 
photoionization. Here are the results. 

1. Photoionization. Putting in the atomic constants and 
calculating the equilibrium temperature and in 
the case of photoionization by a power law of index 
a (F u tx v a ) of either —0.5 or —1, we find that the 
calculated T always lies within the range 0.8-12. The 
explored range in ionization parameter^] U covered all 
the values which produce significant C iv in emission 
(Civ/C > 8%), that is 10~ 3 - 5 < U < 1CT 1 . 

2. Collisional ionization. In this sequence of models, 
we calculated the ionization equilibrium of a plasma 
whose temperature varied from 30 000 K to 50 000 K. 
We find that T remains in the similar low range of 
6-13. At the lower temperature end, Lya emission is 
enhanced considerably by collisional excitation, which 
contributes in reducing T. 

3. Additional heating sources. To cover the case of pho- 
toionization at a higher temperature than the equi- 
librium value (due to additional heating sources such 
as shocks), we artificially increased the photoionized 
plasma temperature to 40 000 K or 50 000 K for calcu- 
lations with the same values of U as above. This did 
not extend the range of T obtained. 

We conclude that for the simple one-zone case, T con- 
sistently remains below the observed value by more than 
two orders of magnitude. 



3.4- The ionization stratified slab 

To verify whether a stratified slab geometry might alter 
the above discrepancy in T, we have calculated in a sim- 
ilar fashion to Bergeron & Stasihska (1986) and Steidel 
(1990b) the internal ionization and temperature structure 
of a slab photoionized by radiation impinging on one-side 
(i.e. one-dimensional "outward only" radiation transfer) 
using the code mappings ic. We adopted a power law 
of index a — — 1 as energy distribution. Since the col- 
umn densities of H and C are useful diagnostics on their 
own right, we present in Fig. g the value of T for a slab 
as a function of Nciv (left panel) and Nhi /Nhii (right 



4 For all practical purposes, the high ionization regime under 
consideration implies that ynu = 1. 

5 We use the customary definition of the ionization parame- 
ter U = ifH/nH as the ratio between the density of ionizing 
photons (impinging on the slab) and the total H density at the 
face of the slab. 
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Fig. 5. a: Calculated and observed T as a function of 
the column density Nciv- The filled circle represents the 
observed value for 0943-242. b: The same models as a 
function of the column ratio Nhi/Nhii- In both pan- 
els, the solid line represents a sequence of photoionized 
slabs with U increasing from left to right, starting at 
10~ 25 . The gas total metallicity is either solar (Z = 1) 
or l/50th solar. The separation between tick marks corre- 
sponds to an increment of 0.25 dex in U. All slab calcula- 
tions were truncated at a depth corresponding to the ob- 
served Nhi = 10 19 cm~ 2 . The slab total column or Nhii 
can be inferred from panel b. [If we were to reduce by 100 
the abundance of the absorption gas while keeping solar 
the emission gas (Z c bs /Z c mi = 0.01, see Eqs | and g), 
this would be equivalent to translating by 2 dex both up 
and to the left the Z = 1 sequence of panel a.] The dotted 
line represents a sequence of slabs of arbitrary uniform 
temperatures (all with U = 10~ 2 and Z = 1) covering 
the range 10 000 K to 40 000 K (from left to right) by in- 
crements of 0.1 dex in T. The open triangle represents a 
slab photoionized by a high velocity shock of V s } loc k = 
500 kms -1 from Dopita & Sutherland 1996. 



panel). (One can interpret Nhi/Nhii of Panel b as the 
mean neutral fraction of the slab: (xjji /uhii) ■) 

The solid line in Fig. || represents a sequence of differ- 
ent slab models with increasing ionization parameter from 
left to right covering the range 10~ 25 < U < 10 _1 for a 
gas of either solar metallicity (Z = 1) or with a signifi- 
cantly reduced metallicity of l/50th solar. The practical 
constraint that C iv be a strong emission line implies that 
U > 10~ 2 ' 5 . In all calculations, the thickness of the slab 
is set by the observable condition that Nhi = 10 19 cm~ 2 . 
Interestingly, such parameters result in a slab which in all 
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cases is "marginally" ionization-bounded with less than 
10% of the ionizing photons not absorbed. 

The monotonic increase of the Nciv column with U is 
in part due to the increasing fraction of C IV but mostly it 
is the result of the slab getting thicker (larger Nhii at con- 
stant Nhi) since xhi decreases monotonically throughout 
the slab with increasing U . The slope or curvature of the 
two solid lines reflect changes in the internal temperature 
stratification of the slab with increasing U. Because of the 
dependence of Y on T (see Eq. ||), there exists an indirect 
dependence of Y on the total metallicity given that the 
equilibrium temperature is governed by collisional excita- 
tion of metal lines (when Z ^> 0.005). 

The striking result from the slab calculations in Fig. || 
is that the models with solar metallicity are still two order 
of magnitudes below the observed Y. Another way of look- 
ing at this discrepancy is to consider separately the 
emission ratio or the ^"iv commn ra tio. Forgetting Y, just 
to achieve the observed column of Nciv (10 145 cm -2 ), 
one would have to use a gas metallicity below solar by 
a factor <; 50 (see sequence with Z = 0.02Zq), which 
cannot be done without irremediably weakening the C iv 
emission line to oblivion. Alternatively, reducing U much 
below 10~ in the solar case can reproduce the Nciv 
column but again the C IV emission line would be totally 
negligible. 

Might the observed t^iv = 0.17 emission line ratio be 
anomalous? This is not the case as the observed value in 
0943-242 is typical of the value observed in others HZRG 
without, for instance, any evidence of dust attenuation of 
Lya. This ratio is also that expected from photoionization 
models if a sufficiently high value of U is used (Villar- 
Martm et al. 1996). 

Another possibility to consider is the presence of other 
heating sources such as shocks which would increase the 
temperature above the equilibrium temperature given by 
photoionization alone. Alternatively, small condensations 
in rapid expansion would result in strong adiabatic cool- 
ing and the temperature would be less than given by cool- 
ing from line emission alone. To explore such cases, we 
have calculated various isothermal photoionized slabs of 
different (but uniform) temperatures (all with U = 10~ 2 ). 
They cover the range 10 000-40 000 K and are represented 
by the dotted line in Fig. [| These models are in no better 
agreement with respect to Y. (Varying U for any of these 
isothermal temperature slabs would result in an horizontal 
line). We also computed Y for a solar metallicity (precur- 
sor) slab submitted to the ionizing flux of a 500 kms -1 
photoionizing shock (Dopita & Sutherland 1996). This 
model which is represented by an open triangle in Fig. ^| 
does not fare better than the power law photoionization 
models. 



4. Discussion 



4-.1. Interpretation of the large Y 

What is the significance of the obvious discrepancy be- 
tween models and the observed Y? Clearly, the working 
hypothesis that the emitting and absorption gas phases 
are physically the same, is now ruled out and an alterna- 
tive explanation must be sought for, based on our result 
that the two gas phases (absorption vs. emission) are phys- 
ically distinct. We consider the two following explanations 
in Sect. [4.1. l| and [4.1. 2| . 



4.1.1. The absorption gas is metal-poor and further out. 

Since the absorption gas in this picture is not spatially 
associated with the emission gas, its metallicity is uncon- 
strained. It turns out that the value of Y ~ 5400 is easily 
reproduced by simply using Z c bs / Z^V 11 ~ 0.005 in the 
one-zone case (see Eqs |J and ^ . The more rigorous strat- 
ified slab geometry would favor a value of Zq/Z^ 1 ~ 
0.01 to reproduce the same T, assuming both gas phases 
to have equal excitation. Can we disentangle the absolute 
abundance values? We cannot rely on the emission spectra 
alone to derive a precise and independent value for Z c bs as 
the emission lines are very model-dependent, with fluxes 
from lines like C iv depending critically on the tempera- 
ture. It can realistically be argued, however, that a Zq" 11 
less than half solar could not reproduce the observed metal 
line ratios. On the other hand, a Z^™ much higher than 
solar cannot be ruled out in absence of direct knowledge 
of the ionizing continuum distribution. We consider more 
plausible a near solar value for Z^V 11 on the ground that 
the extended emission lines extend over 13kpc and there- 
fore sample a huge galactic region very distinct from that 
of the nucelar BLR (hidden here) which has been shown 
to be ultra-solar in high z QSOs (Hamann & Ferland 1999 
and references therein). An attempt, on the other hand, 
to model separately the absorption colu mns observed in 
0943-242 as described below in Sect. 4.2 is more depend- 
able since temperature is much less of an issue. The value 
inferred below of Z c bs ~ O.OIZq is consistent with those 
observed in absorbers of comparable redshift along the 
line of sight of more distant QSOs (Steidel 1990a). Since 
measured galactic metallicity gradients are always nega- 
tive and a function of the distance to the nucleus, such 
a contrast in metallicity between absorption and emission 
gas makes more sense if the absorption gas is located much 
further out than the emission gas which extends to at least 
13kpc in 0943-242. 

We emphasize that this scenario does not entail that 
the absorption gas does not belong to the environment of 
the parent radio galaxy. As argued by v097, the high fre- 
quency of detection of H I aborbers in 9 out of 10 radio 
galaxies smaller than 50 kpc, much in excess of the density 
of absorbers along any line of sight to distant QSOs, is a 
compelling argument for concluding that the absorption 
gas is spatially related to the parent galaxy. Our postu- 
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late is that the large scale H I absorption gas is the same 
gas which is seen instead in emission in those radio galax- 
ies with Lya sizes larger than 50kpc. In effect, absorp- 
tion troughs are not seen when the emission gas extends 
beyond 50kpc. Such objects in general also have much 
larger radio sizes as shown by v097. Kinematically, the 
gas which is seen in emission at the largest spatial scales 
shows narrow FWHM. For instance a reresentative case is 
the radio galaxy 1243+036 (z e = 3.57) which was studied 
in great detail by van Ojik et al. (1996) and which reveals 
the presence of very faint Lya emission extending up to 
136 kpc, a region labelled "outer halo". This emission gas 
has a FWHM of 250 kms -1 and shows clear evidence for 
rotational support. 

A straightforward explanation of why the same gas 
is seen in emission in some objects while in absorption 
in others might simply be the environmental pressure. A 
larger pressure, like the one adopted by v097 can cause 
the warm gas to condense and hence reduce his filling 
factor as compared to similar gas components in a low 
pressure environment. Due to this process, high pressures 
and consequently high densities lead to detectable Lya 
since emissivities scale proportionally to n 2 H: but also to 
an overall smaller covering factor (hence no detectable ab- 
sorption) while low pressures lead to large covering factors 
(hence absorption) as well as negligible emissivities. Dif- 
ferences in pressure in the outer halo would therefore nat- 
urally account for the reported dichotomy of detecting H I 
troughs exclusively in those emission Lya objects devoid 
of very large scale emission (<^ 50 kpc) 

Since absorption troughs tend to be absent in radio 
galaxies showing the largest radio scales, we propose that 
the gas which is seen in absorption must lie outside the 
zone of influence of the radio jet cocoon, a region with 
pressure of order 10 6 Kcm~ 3 (v097). An unpressurized 
outer halo responsible for the absorption troughs ought 
to precede the regime in which the radio material has ex- 
panded sufficiently outward to pressurize the outer halo. 
The eventual increase in environmental pressure would ei- 
ther disrupt the gas or compresses it into small clumps 
(making it unobservable in absorption when the covering 
factor dwindles), which becomes visible in emission if it 
lies within the ionizing cone. v097 assumed that the ab- 
sorption and emission gas were both immersed in zones 
of comparable surrounding pressure (n^T ~ 10 6 Kern" 3 ) 
and were therefore of comparable density (~ 100 cm -3 
for a photoionized gas). We propose instead that when- 
ever aborption troughs are observed, the absorption gas 
must lie outside the radio jet cocoon, allowing for a lower 
density and high covering factor. 

The clear-cut advantages of locating the H I absorber 
in an unpressurized outer halo are threefold: 

1. We can now get the high excitation of the low density 
absorption gas for free. In effect, if the density of the 
absorption gas is as low as 10 _3 -10 -2 cm -3 , the meta- 



galactic background radiation suffices to photoionize 
the absorption gas to the high degree observed in 0943- 
242, whether it does or does not lie within the ionizing 
cone of the nucleus. Conversely, for the objects devoid 
of absorption, when a higher pressure has set in in the 
outer halo (as we presume to be the case in 1243+036), 
the gas is much denser and can be seen in emission only 
if it lies whithin the ionizing cone (since a high density 
gas of ~ 100 cm -3 cannot be kept highly ionized by 
the background metagalactic radiation). This picture 
would be in accord with the findings of van Ojik et al. 
(1996) who detect Lya in emission in 1243+036 only 
along the radio axis (presumably the same axis as that 
of the ionizing radiation cone) and not in the direction 
perpendicular to it. 

2. The much smaller velocity dispersion (b ~ 45 krns" 1 ) 
of the absorption gas as compared to the emission gas 
(FWHM/2.35 ~ 600 kms -1 , cf. Table Q) is more read- 
ily explained if the absorption gas lies undisturbed at 
relatively large distances from the parent galaxy. 

3. It explains why the absorption (yet ionized) gas in 
0943-242 is not seen in emission while being more 
massive than the inner emission Lya gas observed 
within 13 kpc. In effect, the mass of ionized gas ei- 
ther in emission or absorption around 0943-242 in- 
ferred by v097 are 1.410 s M and 10 r (x H i/y H n)- 1 
Mq, respectively. Adopting the conservative value of 
{xhi /vhii) — 0.03 (cf. panel b in Fig. ||), the total 
ionized mass of the absorption ionized gas therefore 
exceed that of the inner emission gas by at least a fac- 
tor two and yet it is not seen in emission! This huge 
pool of ionized gas can remain undetectable in emission 
only if it has a very low density, as argued above. It 
is customary to assume a volume filling factor of 10 -5 
for the gas detected in emission in radio galaxies and 
that this gas is immersed in a region characterized by 
a pressure of order 10 6 K cm~ 3 (v097; van Ojik 1996). 
If we suppose instances where the outer halo has much 
lower pressure than this, it can be shown that for the 
same outer halo mass, the luminosity in Lya would 
scale inversely to the volume filling factor. Hence, the 
gas would be weaker in emission by a factor of 10 if 
its filling factor approached unity (with the mean den- 
sity being lower by the same amount). This scheme 
would easily explain why the outer halo of 0943-242 
is not seen in emission despite its huge mass (com- 
parable incidentally to the outer halo mass measured 
in emission in 1243+036 of 2.8 10 8 Mq by Ojik et al. 
1996). 

4.1.2. A two-phase gas medium 

Due to radiative cooling (which goes as n 2 H and rise 
steeply with T), density enhancements can condense out 
of the emitting gas and form a population of about 100 
times denser and 100 times cooler clouds in pressure equi- 



Binette et al.: A vestige low metallicity gas shell surrounding the radio galaxy 0943-242 at z — 2.92 



9 



librium with the ambient medium. If we maintain that 
the pressure characterizing the absorption and the emis- 
sion gas is comparable (~ 10 6 Kcm -3 ) and that either 
gas phase has a temperature typical of photoionization, 
T ~ 10 4 K, we obtain (adopting a similar notation to v097 
but adapted to the case of 0943-242) that the size and the 
number of small homogeneous absorbing condensations re- 
quired to cover the emission region would be 0.85ro3 pc 
and 2.4 10 8 rQ 3 2 clouds, respectively, where r 3 = 0.03 
x (xhi /vhii)~ 1 [as above we adopt 0.03 as the reference 
neutral H fraction]. Can we find an alternative interpre- 
tation to (1) above for explaining the large V that does 
not require low metallicities for the absorption gas? Such 
a possibility would arise if the Njji column was not di- 
rectly related to the Nciv column. For instance, in the 
auto-gravitating absorber model of Petitjean et al. (1992), 
which consists of a self-gravitating gas condensation with 
a dense neutral core surrounded by photoionized outer lay- 
ers, could in principle give ratios between columns of H I 
and C iv which do not reflect the abundance ratio but rep- 
resents rather the average impact parameter for our line 
of sight. Of course, these models have to be rescaled to 
a pressure of n^T = 10 6 Kcm~ 3 implying much smaller 
sizes but requiring much higher ionizing fluxes (both by 
a factor ~ 10 4 ). This rescaling poses no conceptual prob- 
lems if we assume that the photoionization is by the cen- 
tral AGN. Using their Figures and Table 4 (Petitjean et al. 
1992), we infer that the number of auto-gravitating con- 
densations needed to achieve a covering factor of unity and 
a mean H I column of 10 19 cm -2 would have to be large, in 
excess of 10 95 , for instance, for the model C^og. However, 
after inspection of the various Nciv columns derived from 
their extensive grid of models, we did not find any model 
which would reproduce the observed C IV column without 
having a metallicity < O.IZq. The gain in Z is therefore 
insufficient to get Z^ 

that this explanation for a high T is unworkable 



Z™ ! > 0.5 and we conclude 



4-2. Metallicity determination of the absorption gas 

Our favoured interpretation of the large T is that the ab- 
sorption gas is of very low metallicity compared to the 
(inner/denser) emission gas. Furthermore, a close parallel 
in the physical conditions of the absorption gas could be 
made with those adopted for the study of QSO absorbers 
(e.g. Steidel 1990a, b; Bergeron & Stasinska 1986), namely 
the densities, the metallicities and the excitation mecha- 
nism (photoionization by a hard metagalactic background 
radiation). The observed Njji column of 10 19 cm -2 would 
position the 0943-242 absorber in the category of "Lyman 
limit system" according to Steidel (1992). The coincidence 
in physical conditions might be fortuitous and it does not 
imply per se a common origin or correspondance between 
QSO absorbers and outer halos of radio galaxies. Under 
the sole assumption of similar physical conditions, what 
estimate of the metallicity can we derive for C? From the 



Nciv / Nhi ratio, we cannot determine the ionization pa- 
rameter and therefore directly apply the results and mod- 
els of Steidel (1990b) who determined for each Lyman 
limit system a probable range of U from upper limits or 
from measurements of other species than C iv. It is nev- 
ertheless reasonable to assume that the excitation degree 
in 0943-242 is comparable to that encountered in high 
excitation QSO absorbers. To determine an appropriate 
value for U, we adopted the set of data provided by the 
three Lyman limit systems observed in the spectrum of 
the QSO HS1700+6416 by Vogel & Reimers (1993) who 
successfully measured the columns of up to 3-4 ioniza- 
tion species of each of the three elements C , N and O. 
Amongst our a = — 1 model sequence (Sect. ^4), we se- 
lected the model which had the same U (~ 0.007) as Vogel 
& Reimers (1993) and inferred that the observed columns 
in 0943-242 implied that the Carbon metallicity of the ab- 
sorption gas was 1% solar (that is C/H ~ 410 -6 ), which 
is broadly consistent with the range of Z c bs values favored 
in Sect. |4.1.l| . 



4-- 3. Mean density and cloud sizes 

What would be the minimum density assuming the ab- 
sorption gas to be uniformly distributed? If our pro- 
posed picture was correct, a representative size for the 
absorption gas volume is that given by the outer halo 
as seen in emission in other HZRG. Let us adopt the 
value measured for 1243+036 by van Ojik et al. (1996) 
of 136 kpc. Assuming the same mean ionization parame- 
ter as used above (0.007), we derive a total gas column 
of Nji = Nhi i — 10 21 cm~ 2 . Hence the mean density for 
a volume filling factor unity on a scale of the 1243+036 
outer halo would be ~ 2 10~ 3 cm~ 3 which is a value suffi- 
ciently low to allow photoionization by the feeble ionizing 
metagalactic background radiation. 

4-4- Comparison with the metallicity of BAL QSOs 

Our estimate of the metallicity for the outer halo of 0943- 
242 is at odds with the super-solar metallicities (e.g. 
Hamann 1997, Papovich et al. 2000) of the "associated" 
absorbers seen in high redshift QSOs. The QSO emission 
gas itself (the BLR) is similarly characterized by super- 
solar metallicities (cf. Hamann & Ferland 1999 and refer- 
ences therein). If we consider QSOs and HZRG as equiv- 
alent phenomena observed at different angles, it may ap- 
pear at first surprizing that the metallicities of the ab- 
sorption components are so different. However, we show 
below that this contradiction is only apparent as we are 
probably dealing with totally different gas components. 

1. Kinematics. The HZRG large scale absorbers are kine- 
matically very quiescent. In effect, the modulus of the 
velocity offset between the absorbers and the parent 
galaxy is usually less than 400 km s _1 for the dominant 
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absorber (v097)[J A substantial fraction of HZRG ab- 
sorbers are actually infalling (Binette et al. 1998). 
This is far from being the case for QSO "associated" 
absorbers whose ejection velocities can extend up to 
many thousands kms -1 (Hamann & Ferland 1999). 
For instance, the two associated systems (with de- 
tected metal lines) recently studied by Papovich et al. 
(2000) are blueshifted by 680 and 4900 kms -1 , respec- 
tively. 

2. Selection effect. QSOs are spatially unresolved with 
a size of the source light beam less than a few light- 
weeks across. In the case of HZRG absorbers, the back- 
gound source is the emission gas which extends over a 
scale ~ 35 kpc. This huge difference in scale results in 
a totally different bias on what is preferentially ob- 
served. In effect, the extended absorbers of HZRG are 
weighted towards the largest volumes and hence to- 
wards the most massive gas components (the total 
mass of the absorption component exceeds 10 8 Mq in 
0943-242). By contrast, in the case of QSO associated 
absorbers, the mass of gas directly seen in absorption is 
tiny (e.g. ~ 4 10 -6 Mq if one considers a background 
light beam one light-month diameter and a total gas 
absorption column of 10 18 cm -2 ). 

3. Coexistence with the BLR. To the extent that QSO as- 
sociated absorbers represent gas components expelled 
from the BLR, we should not be surprized that their 
metallicity turn out comparable to the BLR. Given 
that in HZRG we do not directly see the pointlike 
AGN, we cannot expect to see any BLR component in 
absorption. As for the extended gas detected in HZRG, 
there exists no evidence in favour of super-solar metal- 
licities on large scales > 10 kpc (Nv when detected is 
strong only in the nucleus) If a fraction of associated 
absorbers correspond to intervening galaxies close to 
the QSO, we might expect to see amongst counterpart 
HZRGs one or more C iv or Lya absorbers of small 
spatial extent relative to the size of the extended emis- 
sion gas. The weak H I absorption found by Chambers 
et al. (1990) in 4C41.17 might be such occurrence given 
its partial coverage of the Lya background. 

We conclude that HZRG absorbers, when their size is 
comparable to galactic halos (as those found by v097), 
have probably little to do with QSO associated absorbers. 
A more suitable analogy to the absorption gas of HZRG 
is that of the Francis cluster of galaxies at z = 2.38 which 
is characterized by large scale absorption gas on a scale of 
k, 4Mpc (Francis et al. 2000). 

4-5. Constraints on radio galaxy evolution 

The size of the radio source can be used as a clock that 
measures the time elapsed since the start of the radio ac- 

6 Highly blueshifted P-cygni profiles are now known to exist 
in radio galaxies with z > 3.5 (Dey 1999). 



tivity. A number of observed characteristics of distant ra- 
dio galaxies change as a function of radio size, - ie. as 
function of time elapsed (cf. Rottgering et al. 2000). For 
z ~ 1 3CR radio sources, these include optical morphology 
(Best et al. 1996), degree of ionisation, velocity dispersion 
and gas kinematics (Best et al. 2000). At higher redshifts 
(z > 2), only the smaller radio galaxies are affected by Hi 
absorption (v097). All these observations seem to dictate 
an evolutionary scenario in which the radio jet has a dra- 
matic impact on its environment while advancing on its 
way out of the host galaxy (Rottgering et al. 2000, Best 
et al. 2000). 

5. Conclusions 

The detection of Civ absorption in radio galaxy 0943- 
242 at the same redshift as the deep Lya trough observed 
by R095 demonstrates that the detected absorption gas 
is highly ionized. Having assumed that the H I and C iv 
columns measured from the Voigt profile fitting were rep- 
resentative of the dominant gas phase (by mass) in the 
outer halo, we have effectively ruled out that the ab- 
sorption and emission gas occupy the same position in 
0943-242. We subsequently reassessed the picture pro- 
posed by v097 in which both the large scale emission 
gas and the absorption gas were of comparable density 
{riH ~ 100 cm~ 3 ). In the former picture, the absorption 
gas was believed to lie outside the AGN ionization bicone 
(see their Fig. 11 in v097). To ionize the gas to such a 
degree without using the AGN flux is problematic. We 
have proposed an alternative picture in which the absorp- 
tion gas is of very low metallicity and lies far away (in the 
outer halo) from the inner pressurized radio jet cocoon. 
Since in this new scheme the density of the absorption gas 
is expected to be very low, the metagalactic background 
radiation now suffices to photoionizc it. Furthermore, the 
structure of the absorption gas is now drastically simpli- 
fied since we do not need over ~ 10 10 condensations of size 
~ lpc and density ~ 10 2 cm~ 3 to reach a covering factor 
close to unity. We can now reach similarly high covering 
factor using a single or few shells of very low density which 
have a volume filling factor close to unity (assuming a den- 
sity of ~ 10~ 2 ' 5 cm -3 ). 

It appears to us that the low metallicity inferred 
(Z ~ O.OIZq) and the proposed location of the absorp- 
tion gas in 0943-242 -outside the radio cocoon, in an outer 
halo which is seen in emission in other radio galaxies (as 
in 1243+036)- strongly suggest that the absorbers' exis- 
tence precedes the observed AGN phase. Unless this non- 
primordial gas has been enriched by still undetected pop 
III stars, we consider that it more likely corresponds to a 
vestige gas phase expelled from the parent galaxy during 
the initial starburst at the onset of its formation. 

If the C IV doublet was detected in absorption in other 
radio galaxies with deep Lya absorption troughs, there 
are many aspects which would be worth studying. For in- 
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stance, how uniform is the excitation of the absorption 
gas across the region over which it is detected? Is a single 
phase sufficient? This could be tested by an attempt to 
detect absorption troughs of Mgn AA2798 or imaging the 
troughs in C iv with an integral field spectrograph on an 
8-m class telescope. How different is the metallicity of the 
absorption gas in the other radio galaxies? The informa- 
tion gathered could then be used to infer the enrichment 
history of the outer halo gas which surrounds HZRG. 
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